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This paper presents a search for anomalous production of multiple low-energy leptons in associa- 
tion with slWotZ boson using events collected at the CDF experiment corresponding to 5.1 fb~^ 
of integrated luminosity. This search is sensitive to a wide range of topologies with low-momentum 
leptons, including those with the leptons near one another. The observed rates of production of ad- 
ditional electrons and muons are compared with the standard model predictions. No indications of 
phenomena beyond the standard model are found. A 95% confidence level limit is presented on the 
production cross section for a benchmark model of supersymmetric hidden-valley Higgs production. 
Particle identification efficiencies are also provided to enable the calculation of limits on additional 
models. 

PACS numbers: 13.85Qk,12.60Jv,14.80Ly,95.35-l-d 



I. INTRODUCTION 
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The signature of multiple leptons is common in many 
models of physics beyond the standard model (SM) with 
light mass scales and couplings to the electroweak sector, 
such as the Next-to-Minimal Supersymmetric Model [l|, 
little Higgs models 0, and R-parity violating MSSM 
models |3| . Some of these new physics scenarios propose 
explanations for the nature of dark matter Q as well as 
the existence of other, yet-undiscovered particles in long 
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decay chains. In addition to predicting large numbers of 
leptons, these models also often predict that clusters of 
leptons are produced spatially close to each other. These 
clusters are often referred to in the literature as "lep- 
ton jets" 0. Due to the unique characteristics of these 
models, they could have evaded previous searches for an 
excess of leptons, such as diboson searches and SUSY- 
inspircd multi-lcpton searches 0. The high multiplicity 
of leptons can lead to low lepton momenta, well below the 
usual cutoff of 10-20 GeV. Additionally, collimated lep- 
ton jets will fail the standard requirement that leptons be 
isolated in the detector. As an example. Figure [T] shows 
a typical decay chain in a model in which the Higgs de- 
cays to a light hidden sector resulting in events with a 
high multiplicity of leptons [1] . 




luminosity at a center-of-mass energy of ^/s = 1.96 TeV 
collected using the CDF detector at Fermilab between 
December 2004 and January 2010. Within the events 
containing leptonically decaying W and Z bosons, we 
search for additional 'soft' leptons with no isolation re- 
quirements and with momentum _greater than 3 GeV for 
muons and 2 GeV for electrons I 111 . 



II. ANALYSIS STRATEGY 

The analysis strategy and the structure of this paper 
are as follows. The baseline data sets for this analysis 
consist of leptonically decaying W and Z boson events 
selected with high transverse momentum [l^] {pt) lep- 
tons [l^. The kinematic distributions are used to val- 
idate the W and Z boson selections. The selection of 
these events is described in Section IIVI 

After the W ot Z boson reconstruction, additional low- 
Pt electrons and muons are identified in the events with 
no isolation requirements. Purely data-driven techniques 
are used to develop the soft lepton identification algo- 
rithms. The selection of soft leptons is more fully de- 
scribed in Section IVl 

The numbers of additional electrons and muons are 
counted in the inclusive W and Z data sets, where the 
SM predicts few events with multiple leptons. The ob- 
served event count is compared to the SM expectations 
in bins of additional lepton multiplicity. These results 
are described in Section IVlII 



III. THE CDF II DETECTOR 



FIG. 1: An example of multiple low-pr, non-isolated lepton 
production. A Higgs decays to a pair of lightest supersym- 
metric neutralinos (A'^i) which then cascade through a dark 
sector to a lightest dark sector particle (hd) and a number of 
dark photons (7d). The dark photons then decay back into 
the SM in the form of leptons (Z*). This model is adapted 
from Ref. . Note that this diagram shows only the decay of 
the Higgs, while this analysis as a whole would be sensitive to 
the associated production of a Higgs with a, W or Z boson. 



This paper presents a signature-based search for 
anomalous production of multiple electrons and/or 
muons in association with W and Z bosons. Previous 
searches for lepton jets at the Tevatron and at the 
LHC [l3| have focused on searching for clusters of lep- 
tons with specific requirements on the size of the clusters. 
These searches have resulted in no evidence for lepton 
jets. We have performed a more general search, sensitive 
to a wide range of scenarios that predict multiple elec- 
trons and muons. Note that hadronic decays of tau lep- 
tons are not included in this search due to the additional 
difficulty in identifying them in non-isolated topologies. 

The data used here correspond to 5.1 fb~^ of integrated 



The CDF II detector is a cylindrically-symmetric spec- 
trometer designed to study pp collisions at the Fermilab 
Tevatron. The detector has been extensively described in 
detail elsewhere in the literature IJ] . Here the detector 
subsystems relevant for this analysis are described. 

Tracking systems are used to measure the momenta of 
charged particles, to reconstruct primary and secondary 
vertices, and to trigger on and identify leptons with large 
transverse momentum. Silicon strip detectors (SVX) [l5| 
and the central outer tracker (COT) [l^ are contained 
in a superconducting solenoid that generates a magnetic 
field of 1.4 T. The silicon strip system provides up to 8 
measurements in the r — (j) and r — z views and helps to 
reconstruct tracks in the region |r/| < 2 [l2|. The COT 
is an open-cell drift chamber that makes up to 96 mea- 
surements along the track of each charged particle in the 
region \ri\ < 1. Sense wires are arranged in 8 alternating 
axial and ±2° stereo super-layers. The resolution in pT, 
cTp^/pT, is w 0.0015 pT(GeV) for tracks with only COT 
measurements, and « 0.0007 pT(GeV) for tracks with 
both silicon and COT measurements. 

Calorimeters are segmented with towers arranged in 
a projective geometry. Each tower consists of an elec- 
tromagnetic and a hadronic compartment [l7l - [l9| . The 
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central electromagnetic calorimeter (CEM) and central 
hadronic calorimeter (CHA) cover the central region 
(|?7| < 1.1), while the plug electromagnetic calorimeter 
(PEM) and plug hadronic calorimeter (PHA) cover the 
'end plug' region (1.1 < jvyj < 3.6). In this analysis, 
a high-Ex electron is required to be identified in the 
central region, where the CEM has a segmentation of 
15° in (j) and ss 0.1 in rj [l3|, and an Et resolution of 
a{ET)/ET w 13.5%/v/£^T(GeV) © 2% [T^. Two ad- 
ditional systems in the central region with finer spatial 
resolution are used for electron identification. The cen- 
tral strip system (CES) uses a multi-wire proportional 
chamber to make profile measurements of electromag- 
netic showers at a depth of 6 radiation lengths (approx- 
imately shower maximum) [l7| . The central prcshower 
detector (CPR) is located just outside the solenoid coil 
on the front face of the CEM. In 2004 the CPR was up- 
graded from the Run I configuration of wire proportional 
chambers to a fast scintillator system [l^. This analysis 
only uses data collected after the CPR upgrade. 

Muons arc identified using the central muon sys- 
tems [23 : CMU and CMP for the pseudo-rapidity region 
of I77I < 0.6, and CMX for the pseudo-rapidity region of 
0.6 < I77I < 1.0. The CMU system uses four layers of pla- 
nar drift chambers to detect muons with > 1-4 GeV. 
The CMP system consists of an additional four layers of 
planar drift chambers located behind 0.6 m of steel out- 
side the magnetic return yoke, and detects muons with 
PT > 2.2 GcV. The CMX system detects muons with 
PT > 1.4 GeV with four to eight layers of drift chambers, 
depending on the direction of the muon. 

The luminosity is measured using two sets of gas 
Cerenkov counters [2l[, located in the region 3.7 < \ri\ < 
4.7. The total uncertainty on the luminosity is estimated 
to be 5.9%, where 4.4% comes from the acceptance and 
operation of the luminosity monitor and 4.0% from the 
calculation of the inelastic pp cross-section [22| ■ 

A three-level online event selection (trigger) sys- 
tem [2^ selects events to be recorded for further anal- 
ysis. The first two trigger levels consist of dedicated fast 
digital electronics analyzing a subset of the complete de- 
tector information. The third level, applied to the full 
set of detector information from those events passing the 
first two levels, consists of a farm of computers that re- 
construct the data and apply selection criteria consistent 
with the subsequent offline event processing. 



IV. W AND Z BOSON SAMPLE SELECTION 

Events for this analysis are selected with three different 
triggers [23|. Approximately half the events are selected 
with a trigger requiring a high-p^ central electron in the 
CEM {Et > 18 GeV, |?7| < 1.0). In addition, two muon 
triggers, one requiring hits in both the CMP and CMU 
and the other requiring hits in the CMX, collect events 
with central muons {px > 18 GeV, |?7| < 1.0). 

Further selection criteria are imposed on triggered 



events offline. Electron (muon) candidates are required 
to have Et > 20 GeV {pr > 20 GeV). They must ful- 
fill several other identification criteria designed to select 
pure samples of high-px electrons (muons) |13| , including 
an isolation requirement that the energy within a cone of 
AR = ^ A(/)2 + A772 < 0.4 around the lepton direction is 
less than 10% of the Et [pt) of the electron (muon). 

In order to reduce the electron background from pho- 
ton conversions, the electron(s) from the W or Z boson 
decay are required to pass a conversion filter. Electron 
candidates with an oppositely-charged partner track con- 
sistent with having originated from a photon conversion 
are removed |2J|. However, the electron candidate is kept 
if its partner conversion track also has another partner 
track, since the three tracks are assumed to originate 
from an electron which radiates a photon which subse- 
quently converts. 

In order to reduce the background from mesons decay- 
ing to muons within the tracking chamber, the muon(s) 
from the W 01 Z boson decay must pass a decay-in-fiight 
(DIE) removal algorithm. The DIE algorithm requires 
the P^r degree of freedom of the fitted track to be 
less than 3 and requires that the impact parameter of 
the track be less than 0.02 cm. Additionally, for tracks 
with Pt > 300 GeV, it requires A'transitions > 30, where 
-^transitions is the number of times the pattern of track 
hits crosses the fitted track [2^. Muons consistent with 
cosmic rays are vetoed [2^ . 

To select W boson events we require ^t > 25 GeV 
and that the highest-energy lepton and the ^t have 
niT > 20 GeV In order to remove events where 

the ^T arises from a mismeasured lepton, the difference 
in (j> between the highest-energy lepton and the ^t is 
required to be greater than 0.5 radians. The Z boson se- 
lection requires two oppositely-charged, same-flavor lep- 
tons. One of these leptons is required to pass the above 
high-pT lepton identification selections while the other is 
required only to pass a less stringent 'loose' selection. For 
muons. the loose selection allows for muons with pT > 10 
GeV that have hits in either the CMP, CMU, or CMX 
systems. For electrons, the loose selection accepts elec- 
trons with Et > 12 GeV and has relaxed identification 
requirements with respect to the centroid shape in the 
CES and E/p, the ratio of calorimeter energy to track 
momentum |27| . Finally, the invariant mass of the lep- 
ton pair is required to be within the range of 76 GeV 
< m{l,l) < 106 GeV, consistent with the mass of the Z 
boson. 

The distributions of my in W boson events and the 
dilepton invariant mass in Z boson events arc shown in 
Figure [2] for both electron- and muon-triggered events. 
In total, 4,722,370 W boson events and 342,291 Z bo- 
son events are obtained from 5.1 fb~^ of data. Good 
agreement with predictions is observed across most of 
the distributions. In the W uit distributions, a disagree- 
ment occurs at low mass, where the distribution shifts 
from being QCD-dominated to electrowcak-dominated, 
and is accounted for by the QCD normalization system- 
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atic uncertainty (as described in Section I VI A[) . In the 
Z selection, a similar mass disagreement is due to the 
fact that the Monte Carlo (MC) simulation does not in- 
clude Drell-Yan events with a Z j^* invariant mass below 
8 GeV. It is eliminated with the requirement that the 
dilepton mass be within the Z peak. 



V. SOFT LEPTON IDENTIFICATION 

The identification of low-pr, or "soft", Icptons is a 
main focus of this analysis. Likelihood-based methods 
are used to identify soft electrons and muons. The iden- 
tification algorithms are described here, along with the 
methods used to validate them and evaluate their sys- 
tematic uncertainties. 



A. Soft electrons 

Soft electrons arc identified using a likelihood method 
trained on a signal sample from photon conversions and 
a background sample from other tracks with electron 
sources removed. 



1. Identification algorithm and candidate selections 

A preselection is applied to all soft electron candidates 
requiring good track quality as well as track extrapola- 
tion to the CES, CPR, and calorimeter. Only tracks with 
jryl < 1 are considered for the soft electron identification. 

After this preselection, a likelihood-based calculator 
is used to identify electrons. The likelihood calculator 
uses seven discriminating variables: the energy loss as the 
track traverses the tracking chamber, the electromagnetic 
and hadronic calorimeter energies, the energies deposited 
in the preradiator and the showermax detector, and the 
two-dimensional distance (Ax, Az) between the extrap- 
olated position of the track and the shower in the CES. 
The calorimeter variables are calculated using a narrow, 
two-tower- wide section of the calorimeter. 

Some of the variables used in the soft electron identifi- 
cation are modeled very badly in the MC, and so the like- 
lihood is trained on data without resorting to the simula- 
tion. For each of the above variables Xi, a fit is performed 
to the ratio of the distribution in the electron sample and 
the distribution in the non-electron background ( "fake" ) 
sample. For each candidate, the value of each of these fit 
functions is multiplied together to get the final likelihood 

(-^electron) ■ 
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The distribution of the likelihood in the real and fake 
samples is shown in Figure [31 A candidate is identified as 
an electron if it passes the requirement ^electron > 0.99. 
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FIG. 2: (a) The transverse mass (mr) of the highest-pT lepton 
and the in the electron-triggered W boson sample, (b) 
The dilepton invariant mass in the electron-triggered Z boson 
sample, (c) The mr of the highest-pT lepton and the 
in the muon-triggered W boson sample, (d) The dilepton 
invariant mass in the muon-triggered Z boson sample. The 
estimation of the QCD contribution to these distributions is 
described in Section TVI Al The points represent the observed 
data and the filled histograms are the SM estimates. 
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FIG. 3: (a) The likelihood distributions for electrons (closed 
squares) and non-electrons (open circles) after all preselection 
criteria. Only those candidates with a likelihood > 0.99 are 
identified as electrons, (b) The efficiency as a function of pr 
for the identification of electrons (closed squares) and tracks 
misidentified as electrons (open circles) after the likelihood 
selection. 



Training samples and efficiency and misidentification 
rate measurements 



Photon conversions are used as a pure sample of elec- 
trons to train the likelihood function. In events selected 
using an 8 GeV electron trigger, pairs of tracks are found 
that correspond to a photon converting into e+e~ [23 |. 
In order to avoid any bias from the trigger, the lower- 
momentum track of the conversion pair is used to train 
the likelihood. 

Events from the 18 GcV muon trigger arc used to select 
a sample of non-electron tracks with which to train the 
likelihood function. All tracks in the events that, along 
with another track, form a possible photon conversion are 
removed from the training sample. To reduce the bias 
from using a muon-triggered sample, any track that is 
within AR < 0.7 of an identified muon is also ignored. In 
addition, to reduce the contamination from real electrons, 
any event that contains an identified heavy quark decay 
or an identified high-pT electron is ignored. 

The efficiency and fake rate are calculated in these 
training samples as functions of pT, i], and track iso- 



lation. The same sample used for training is also used 
to measure the efficiency, due to the larger backgrounds 
present in other independent samples. The separation in 
identification rate between electrons and non-electrons 
after the likelihood selection is shown in Figure [3] (right). 
The efficiency in terms of px and 77, after the track and 
CES shower have been identified, is shown in Table HI 



TABLE I: Efficiency to identify soft (2 GeV < pt < 20 GeV) 
electrons as a function of candidate pr and 77. 

PT range (GeV) [2, 2.5] [2.5, 3] [3, 6] [6, 12] [12, 20] 

< |r?| < 0.2 0.90 0.99 0.99 0.90 0.99 

0.2 < \ri\ < 0.6 0.94 0.99 0.99 0.99 0.99 

0.6 < |r?| < 1 0.95 0.99 0.99 0.99 0.99 



This identification rate is applied as a weight to each 
candidate track in the MC to find the predicted number 
of identified electrons. 



3. Validation and systematic uncertainty determination 

The efficiency and fake rate parametrizations are 
checked on a data set triggered on jets having Et > 50 
GeV. The parametrizations use the px, Tj and isolation 
of candidates in order to account for any kinematic dif- 
ferences between the training sample and the validation 
sample. First, the same electron removal that was used 
for the fake training sample (Section IV A2[) is applied to 
the tracks in the jet sample. The likelihood distribution 
of all candidate tracks in the jet sample is then fit to tem- 
plates from the real and fake likelihood training samples 
to obtain the fraction of real and fake electrons in the 
jet sample. The jet sample is found to consist of 2.5% 
real electrons, mostly coming from photon conversions 
from which only one electron was reconstructed. The 
predicted identification rate is then checked for agree- 
ment with the measured identification rate. 

The disagreement between the calculated and observed 
identification rates is measured to be 1.6%. However, we 
observe larger disagreement in the shapes of the calcu- 
lated and observed distributions in pT and 77. We assign a 
systematic uncertainty of 15%, which is sufficient to cover 
the observed disagreement [2J]. This systematic uncer- 
tainty is applied separately to the electron identification 
and misidentification rates. 



B. Soft muons 

Soft muons are identified using a method similar to 
that described in Ref. [2^. The inputs to the algorithm 
are derived from a sample of muons arising from J /ip 
decays in a muon calibration dataset. 
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1. Identification algorithm and candidate selections 

The soft muon identification algorithm relies on match- 
ing tracks identified in the COT to track segments recon- 
structed in the muon chambers (muon stubs). Match- 
ing is done in the extrapolated position along the muon 
chamber drift direction {x), the longitudinal coordinate 
along the chamber wires {z) when available, and the dif- 
ference in slope between the extrapolated COT track and 
the reconstructed muon chamber track segment {4>l)- 
Tracks are paired with muon chamber track segments 
based on the best match in x for those track segments 
within 50 cm of an extrapolated COT track. 

Soft muon candidates arc required to extrapolate to 
within the physical boundaries of a muon chamber, have 
good track quality, have at least one hit in the SVX, 
|c?o| < 0-3 cm where do is the impact parameter with 
respect to the beamline, and zq < 60 cm where zg is the 
z position of the track at the interaction point. 

A is built from the track-to-stub matching variables 
Xi described above (da;, dz, and dfj)^). This is normal- 
ized to have mean and variance 1 for real muons, inde- 
pendent of the number n of track-stub matching variables 
Xi used: 



E 



(Xi - ^.lif 



where fii and erf are the expected mean and variance of 
the distribution of x^, and (t^x^) is the expected standard 
deviation of x^- 

In the final selection, we require that all identified soft 
muons must have a track segment in each muon chamber 
to which the track extrapolates and that |Qmuon| < 3.5 
(see Fig. SI). 



2. Efficiency and misidentification rate measurements 

The efficiency of the soft muon identification is mea- 
sured using a pure sample of muons obtained from J/ijj — > 
/i/j, decays. These events are obtained using an online 
trigger requiring the presence of a muon with pT > 8 
GcV. The J/tp is reconstructed by requiring that the 
trigger muon make a vertex with another track of op- 
posite charge that has associated muon chamber hits. 
All track requirements listed in Sec. IV B Tj are applied to 
both tracks. The J/'^ candidate mass is required to sat- 
isfy 3.03 < m{fifj,) < 3.15 GeV, and the sidebands of the 
mass distribution are used to evaluate the background 
under the mass peak. 

The misidentification rates of pions and kaons are mea- 
sured in Z3*+ £'°7r+ decays where the decays 
as — > K^TT^. These events are obtained from a 
trigger that requires the presence of a vertex contain- 
ing two tracks and are reconstructed requiring masses 
1.835 < m{KTT) < 1.895 GcV and m{D*) - m{D°) < 170 



McV. The sidebands of the m{D*) — m(D°) distribution 
are used to evaluate the background under the mass peak. 

The misidentification rate of protons is measured us- 
ing a sample of protons obtained from A — p7r decays. 
These events arc taken from the same dataset as that 
from which the D* sample is obtained. The reconstructed 
A mass is required to satisfy 1.111 < m{pn) < 1.121 
GeV. The sidebands of the mass distribution are used to 
evaluate the background under the mass peak. 

Figure S] (left) shows the distribution of muon scaled 
1, using the samples described above. Good 
separation is obtained between muons and other particle 
species. 
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FIG. 4: (a) The distribution of soft muon scaled x^i Qmuon, 
for muons, pions, kaons, and protons after all preselection se- 
lections. Only those candidates with | Qmuon | < 3.5 are identi- 
fied as muons. (b) The muon identification rate (circles) and 
misidentification rates for pions (triangles), kaons (squares), 
and protons (triangles) after the scaled selection. 



An efficiency matrix is created in bins of pT and rj 
using the J/^j sample. Because the sample is limited 
in statistics for px > 12 GeV, empty bins are filled in 
using interpolation between the low-px muons from J/^P 
decays and higher-pr muons from Z decays. The soft 
muon identification is applied to Z events so that the 
region between the J/ip and Z px may be correctly fitted. 
Note that Figure |3] shows the observed results in these 
low-statistics bins, while Table HIl shows the interpolated 
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TABLE II: Efficiency to identify soft muons as a function of candidate pr and 77. 



PT range (GeV) [3, 4] 

-1.5 < 77 < -0.7 0.739 

-0.7 < 77 < -0.55 0.593 

-0.55 < 77 < -0.45 0.749 

-0.45 < 77 < -0.15 0.816 

-0.15 < 77 < 0.15 0.777 

0.15 < 77 < 0.45 0.832 

0.45 < 77 < 0.55 0.768 

0.55 < 77 < 0.7 0.625 

0.7 < 77 < 1.5 0.750 



[4,5] [5,6] [6,8] [8,10] 

0.626 0.567 0.419 0.342 

0.556 0.581 0.480 0.438 

0.788 0.883 0.751 0.783 

0.901 0.898 0.782 0.821 

0.796 0.784 0.667 0.657 

0.918 0.913 0.799 0.815 

0.782 0.840 0.741 0.582 

0.573 0.556 0.461 0.450 

0.617 0.593 0.428 0.327 



[10,12] [12,16] [16,20] 

0.127 0.237 0.174 

0.299 0.356 0.344 

0.608 0.644 0.659 

0.701 0.570 0.659 

0.525 0.424 0.616 

0.698 0.568 0.659 

0.758 0.529 0.659 

0.409 0.237 0.256 

0.146 0.173 0.174 



efficiencies. 

For the corresponding binned misidentification matrix, 
the misidentification rate is measured in each of the three 
background samples. The n, K, and p matrices are then 
combined in the proportion found in W boson decays. 
These relative proportions arc found to be /(tt) = 0.719, 
f{K) = 0.156, and f{p) = 0.125. 

The efficiency and fake rate as a function of pT is shown 
in Figure |31 The efficiency in terms of pT and 7/ is tabu- 
lated in Table HIl These identification rates are measured 
after the track and muon hit(s) have been identified. 

This identification rate is applied as a weight to each 
candidate track in the MC to find the predicted number 
of identified electrons. 



3. Soft muon systematic uncertainty determination 

Separate systematic uncertainties are estimated for the 
true muon identification efficiency and the misidentifica- 
tion rate. The invariant mass sideband subtraction tech- 
nique used to obtain the muon efficiency matrix intro- 
duces uncertainties arising from the statistics of the J/ 
sample. These uncertainties vary from 2% - 70%, depend- 
ing on the bin in pT and 77. In addition, the maximum 
variation in efficiency of 8% arising from the difference 
between isolated and non-isolated candidates is used as 
an uncertainty representing the maximum possible differ- 
ence between the J/ ip sample environment and the W/ Z 
environment. This is added in quadrature to the statis- 
tical uncertainty arising from the sideband subtraction 
method to obtain a final muon efficiency uncertainty of 
8% - 70%. 

The misidentification systematic uncertainty is ob- 
tained by selecting muon-free regions in samples triggered 
on high-pr jets and taking the difference between ob- 
served and predicted soft muon misidentification rates. 
In this jet sample, at least 3 jets are required with 
Et > 15 GeV and 177] < 2.0. In order to reduce the 
contamination from real muons, any jet that contains an 
identified heavy quark decay is rejected, as is any track 
that has impact parameter significance do/a{do) > 2. In 
a sample having an online trigger requiring the presence 



of a jet with Et > 100 GeV, a difference of 4.8% is 
observed between the observed and predicted soft muon 
identification rates. A conservative estimate of twice this 
difference is used as the systematic uncertainty on the 
soft muon misidentification rate. 

C. Application of soft lepton identification to W/Z 
samples 

Additional selection criteria are applied to soft lepton 
candidates in the high pT W and Z boson data samples 
to reduce the amount of background in the search sam- 
ple. Any track that is already identified as a high-pr 
electron or muon in the W or Z boson selection is inel- 
igible to be identified as a soft muon. To reject badly 
measured tracks, each track is required to have at least 
one hit in the silicon detector. For electron candidates, 
this hit is required to be within the first two layers of the 
silicon detector to help reject photon conversions. Each 
track is required to be inside of a reconstructed jet hav- 
ing ]77] < 2.0 and transverse energy of Et > 5 GeV, so 
that the heavy fiavor fraction fit described later in Sec- 
tion IVIBI can be applied. (Note that the 'jet' could be 
composed entirely of leptons, or even entirely of a sin- 
gle lepton.) Any track that is identified as a conversion 
partner is rejected. The track candidate must have a dis- 
tance along the beamline \A.z\ < 5 cm from the high-pr 
trigger lepton. If the trigger lepton is the same flavor as 
the soft lepton, the invariant mass AI is calculated of the 
candidate and trigger, and the following mass ranges are 
rejected: 

• A/ < 5 GcV to suppress the J/ip and bb back- 
grounds. 

• 9 < A/ < 10 GcV if the candidate track has op- 
posite charge to the trigger lepton. This rejects T 
events. 

• 80 < M < 100 GeV if the candidate track has 
opposite charge to the trigger lepton. This rejects 
Z events. 

These additional selection criteria have a small effect 
on the benchmark model chosen for this analysis, cutting 
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out 4.5% of the signal leptons generated. 



VI. BACKGROUND PREDICTION 

The main SM backgrounds in this analysis are from 
W + jets, Drell-Yan, QCD multijet, top quark, and di- 
boson production processes. The cross section and dif- 
ferential distributions of electroweak backgrounds from 
hard scattering processes are modeled using the Alp- 
GEN [2§| MC program, except for the top production and 
diboson production backgrounds, which are modeled by 
Pythia [301 ■ Pythia is used to model the parton show- 
ering in all samples. These MC events are analyzed using 
a GEANT based detector simulation [3l|. The samples 
generated by Alpgen are W/Z + Np partons (light fla- 
vor) and W + qq + Np partons, where q = c,b (heavy 
flavor). The interface with the parton showering gener- 
ates a double counting of heavy flavor events, which is 
corrected using the MLM matching method [S^] ■ 

The relative contributions from the various back- 
ground sources can be seen qualitatively in Figured The 
cross sections used for every sample are described in . 
The final background predictions are summarized later 
in Section IVlII The QCD multijet background requires a 
different treatment since it is not possible to simulate it 
using MC. It is derived using data as explained below. 



A. QCD multijet background fraction 



GeV. The MC electroweak contribution and the data- 
derived QCD template are scaled to the result obtained 
from this fit. Figure[S]shows the result of this fit in the 
electron-triggered dataset. A similar fit is performed in 
each muon-triggered dataset. A systematic uncertainty 
of 26% is applied to the QCD normalization, as found 
in M. 
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FIG. 5: The fit to the distribution of events with mr > 20 
GeV and A(^(_^t,0 > in the electron-triggered dataset. 
The "electroweak" template is obtained from Monte Carlo 
and the "QCD" template is obtained from the anti-selected 
electron data sample. The systematic uncertainty of 26% 
found in [331 is shown. 



The W boson is identified by the presence of a high 
energy lepton and missing transverse energy. Events con- 
taining jets may emulate this signature; a dijet event, for 
example, may have large arising from the energy mis- 
measurement of one jet while the other jet in the event 
can mimic an electron by leaving a track in the COT 
associated with an electromagnetic energy deposit. The 
contribution from these QCD multijet processes is esti- 
mated by using a data-derived model [331 ■ This is ac- 
complished by defining an object that is similar to an 
electron, but has a much larger rate of contamination 
from jets; we refer to this as an "anti-selected electron" . 
An anti-selected electron is required to pass the same 
kinematic requirements as an electron, but must fail at 
least two of the identification requirements. 

The number of events arising from the QCD multijet 
background is obtained by fitting the distribution of 
the data using two templates: an electroweak template 
obtained from W+ jets, Z+ jets and diboson MC, and 
a QCD template. The QCD template is obtained from 
the anti-selected electron sample after subtracting the ex- 
pected W boson contamination using the MC. The total 
number of events is kept constant and the fraction from 
each template is obtained from the fit. 

After the fit is performed across the distribution, 
the number of QCD events in the W boson signal re- 
gion is calculated by applying the selection of > 25 



B. Heavy flavor background fraction 

The leptonic decay of heavy fiavor quarks creates a 
significant background contribution to the soft leptons 
of this analysis. This background is estimated using the 
data in the W/Z + exactly one soft muon channel, which 
should be dominated by SM processes. A fit is performed 
in two distributions of soft muons which are sensitive to 
the heavy flavor fraction: p^^' , which is the momentum of 
the muon transverse to the direction of the jet in which 
it is found, and do/a{dQ), which is the significance of 
the muon's impact parameter with respect to the beam- 
line. A simultaneous fit is performed of these two dis- 
tributions to a sum of templates from heavy fiavor, light 
flavor, and Drell-Yan processes, as shown in Figure |6l 
These templates were acquired from the MC background 
samples. The result of this fit is used to normalize the 
contributions of the three types of processes in the higher- 
multiplicity sample. The uncertainty resulting from the 
fit. ranging from 5% to 34% in the various samples, is 
used as a systematic uncertainty on this normalization. 
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FIG. 6: The result of the simultaneous fit of the + 1 soft muon sample in the and do significance of the soft muon. The 
data distribution is fit to the sum of three components: PF+heavy quark, W+light quark/gluon, and Drell-Yan. 



C. Normalization of soft electron multiplicities 

The heavy flavor fit described in Section IVIBI normal- 
izes all of the data to the W/Z+lfj, channel. However, 
we find a mismatch in the W / Z+le channel, which has 
a large contribution from photon conversions. The dif- 
ference between the predicted and observed numbers in 
the W/Z plus exactly one electron channel is 34% in the 
W boson sample and 31% in the Z boson sample. This 
is used as a systematic uncertainty for the normalization 
of all other MC with at least one additional identified 
electron (23 |. 

VII. RESULTS 

Using the soft lepton identification techniques de- 
scribed in Section |Vl we count the numbers of W and 
Z boson events with multiple additional leptons. Fig- 
ures [7] and [5] show the multiplicity of additional electrons 
(iVe) and muons (N^) in these events, with the SM expec- 
tation and observed data overlaid. The two-dimensional 
histograms of Nfj^ vs. are presented in slices of 
for ease of viewing. These expected and observed event 



counts are also presented in Tables IIIII and IIVI for ease 
of comparison with predictions from other models. The 
sources of systematic uncertainties are summarized in Ta- 
ble |Vl with references to the sections in which they are 
described and evaluated. Good agreement with the SM 
expectation is observed across the distributions. 

In particular, very few multi-muon events are observed. 
This is the region where many lepton jet models would 
be expected to show an excess, since a potential signal 
in the multi-electron region would be more likely to be 
hidden by the large background contribution from pho- 
ton conversions. Only three events containing 3 muons 
beyond the W selection are observed, which is consistent 
with the SM expectation of 2.9 events. No events are 
observed containing four or more additional muons. 



A. Benchmark model 

This is a general signature-based search, and as such 
is applicable to many different models. We choose an ex- 
ample model from the representative lepton jet models 
presented in Ref. The benchmark model chosen for 
this analysis is an adaptation of the 'Neutralino Bench- 
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FIG. 7: Multiplicity of additional electrons and muons after the W boson selection. The two-dimensional histogram of A^^ vs. 
Afg is presented in slices of A'^e for ease of viewing. Both hard and soft leptons (but not the initial lepton used for the W boson 
selection) are counted. Note that the distributions combine the electron- and muon-triggered events. 



mark Model,' in which the Higgs decays principally to a 
pair of the lightest supersymmetric particles, which then 
decay through a dark sector to lepton jets. A MC sam- 
ple of signal events was generated from this model using 
Pythia. The signal from this model to which this anal- 
ysis is most sensitive is associated production of a or 
Z boson and a Higgs boson, which has a cross section of 
389 fb. This cross section would result in 1647 Higgs 
events and 322 Higgs events in the data sample of 



this analysis before applying any selection criteria. 



The particular parameters of the model [3J] were cho- 
sen to create a 'typical' model of this class. The MSSM 
parameters (^, mi, m2, tan(/3) and sin(a)) avoid previ- 
ous limits from searches for supersymmetry while mak- 
ing the lightest supersymmetric partner (xo) the favored 
Higgs decay channel. The Higgs has a mass near that 
favored by precision measurements. The branching frac- 
tions for xo decaying into the dark ncutralinos (x^) and 
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FIG. 8: Multiplicity of additional electrons and muons after the Z selection. The two-dimensional histogram of A'^ vs. Ne 
is presented in slices of for ease of viewing. Both hard and soft leptons (but not the initial leptons used for the Z boson 
selection) are counted. Note that the distributions combine the electron- and muon-triggered events. 



dark photons (7^) simply model the sort of cascade decay 
illustrated in Figure [T] The mass of the dark photon is 
chosen in order to make the additional leptons that are 
produced approximately half muons and half electrons. 
These parameters are summarized in Table IVIl 

We set a 95% confidence level limit on the production 
of this benchmark model. The limit is set at 0.312 x ct, or 
112 fb. The model can be ruled out at the standard cross 
section at a confidence level of 99.7%. Both of these limits 



are set in the Bayesian framework using the mclimit 
tools [35! running over the combined W and Z channels 
in Figures [7] and [5] (Tables lllll and HV)) . 

B. Application to other models 

In addition to the benchmark model discussed in Sec- 
tion IVII Al limits can be set on a wide range of alternate 
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TABLE III: Summary of predicted and observed event counts by number of additional electrons (A*'e) and muons (Nfj,) after 
the W boson selection. The prediction of a model described in Section FVII Al ls also shown for comparison. Bins with less than 
0.25 expected events in both signal and background and observed events are not shown. 



Ne Predicted SM Background Predicted Dark Higgs Signal Observed 
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0.015 






models. A rough estimate of the limit for a particular 
model can be made by normalizing its production to the 
W 01 Z boson cross section, applying the efficiencies in 
Tables H] and |lT] to the additional leptons, and compar- 
ing the result to the observed and predicted numbers of 
additional leptons in Tables IIIII and IIVI For ease of refer- 
ence, a summary of the kinematic selections for identified 
objects is presented in Table IVlIl 



VIII. CONCLUSIONS 

This analysis expands the reach of previous searches 
for additional leptons by allowing leptons to be recon- 
structed from a much lower pT threshold and with no 
requirement of isolation. This greatly increases the ac- 
ceptance to find Icpton jets or similar excesses of leptons 
from effects beyond the SM. No indication of such new 
effects is seen in the data sample. A 95% confidence level 
limit is set on an example benchmark model of supersym- 
metric Higgs production, and a framework is provided to 
set limits on a class of other models. 



In general, any model that predicts significant numbers 
of 3-muon events can be ruled out, since only three such 
events are observed in the sample, consistent with the 
SM background. However, models that produce multiple 
electrons can more easily be accommodated, since photon 
conversions result in a much higher background in that 
region. 
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TABLE IV: Summary of predicted and observed event counts by number of additional electrons (Ne) and muons (Nf^) after 
the Z selection. The prediction of a model described in Section [VII Al is also shown for comparison. Bins with less than 0.25 
expected events in both signal and background and observed events are not shown. 



Ne Predicted SM Background Predicted Dark Higgs Signal Observed 



1 
1 1 
1 2 



2 

2 1 

3 

3 1 

4 

5 



215219 ± 36886 
255 ± 52 
3.2 ±0.89 

2145 ± 447 

30 ±8.1 
0.51 ±0.18 

175 ± 50 
4.2 ± 1.5 
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0.52 ±0.35 
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0.070 
0.031 

0.019 

0.0044 



211448 
270 
4 

1975 
20 
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TABLE V: Sources of systematic uncertainties. Their size is measured both as a percentage and as the number of events in a 
benchmark-signal-rich region, defined as a or ^ boson plus at least 3 additional muons withpT > 3 GeV. Note that, although 
some of the systematics are large, they have little effect in the signal region due to there being negligible SM background. 



Systematic Source 



Size 



Effect in Large S/B Region (Events) 



Trigger Efficiency [ 24] ±(1.6 - 5.9)% 

QCD fraction (Sec. WT^i ±26% 

Soft e real rate (Sec. IV A 31) ±15% 

Soft e fake rate (Sec. IVX3ll ±15% 

Soft real rate (Sec. IVB3|l ±(8- 70)% 

Soft fi fake rate ('Sec. lVB3|l ±10% 

Soft e normahzation (Se dVICl) ±(31-39)% 

Heavy Flavor Fraction (Sec. NlB^ ±(5-34)% 



±0.06 


±0.04 
±0.11 
±0.64 
±0.34 
±0.24 
±0.25 



TABLE VI: Parameters used for the benchmark model based 
on that in Ref. [8[. The first five parameters are the inputs to 
the MSSM including the branching fractions for Xo Xd + 



Parameter 



Value 



mi (bino) 
7712 (wino) 
tan(/3) 
sin(a) 

ruH 



m 



'"Id 



BR(xo ^Xd + 27^; 
BR(xo -s- Xd ± 37^; 
BR(XQ ^ Xd ± ^Id] 



149 GeV 
13 GeV 
286 GeV 
3.5 

-0.28 
10 GeV 
120 GeV 

1 GeV 
300 MeV 

33% 

33% 

33% 
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TABLE VII: Summary of kinematic requirements to find various objects. These numbers can be used to set limits on many 
models that predict production of additional leptons. 



Object Requirements Number Observed 

PT(e/^i) > 20 GeV 4,722,370 
|»7(e)| < 1.1, < 1-5 

> 25 GeV 
mT{l,:pT) > 20 GeV 
d(j>{l,^T) > 0.5 

Z PT{e/fi) > 20 GeV 342,291 

PT(e2) > 12 GeV, pt(m2) > 10 GeV 
\v{e)\ < 1.1, < 1.5 

76 GeV < m(li,l2) < 106 GeV 

soft e pT(e)>2GeV See Tables IIIIl and HVl 

\v{e)\ < 1 
£ > 0.99 (Efficiency in Table IJ) 

soft pT(/i)>3GeV See Tables Hn] and [IV] 

< 1-5 

|£| < 3.5 (Efficiency in Table [Il| 



[1] R. Dermisek and J. F. Gunion, Phys. Rev. Lett. 95, 
041801 (2005); R. Dermisek and J. F. Gunion, Phys. Rev. 
D. 75, 075019 (2007). 

[2] C. T. Hill and R. J. HiU, Phys. Rev. D 75, 115009 (2007); 
Phys. Rev. D 76, 115014 (2007). 

[3] C. Cheung, J. T. Ruderman, L. Wang, and I. Yavin, J. 
High Energy Phys 04, 116 (2010). 

[4] N. Arkani-Hamed, D. P. Finkbeiner, T. R. Slayter 
and N. Weiner, Phys. Rev. D 79, 015014 (2009); 
P. Meade, M. Papucci, A. Strumia, and T. Volansky, 
arXiv:0905.0480 [hep-ph]; M. Pospelov and A. Ritz, 
Phys. Lett. B 671, 391 (2009); C. Cheung, J. T. Ruder- 
man, L. T. Wang and I. Yavin, Phys. Rev. D 80, 035008 
(2009); A. Katz and R. Sundrum, J. High Energy Phys. 
06, 003 (2009). 

[5] J. T. Ruderman and T. Volansky, arXiv:0908.1570 [hep- 
ph]; X. Chen, JCAP 0909, 029 (2009); J. Mardon, Y. No- 
mura and J. Thaler, Phys. Rev. D 80, 035013 (2009). 

[6] T. Aahonen et al. (CDF Collaboration), Phys. Rev. D 
83, 112008 (2011). 

[7] T. Aaltonen et al. (CDF Collaboration), Phys. Rev. Lett. 
101, 251801 (2008); V.M. Abazov et al. (DO Collabora- 
tion), Phys. Lett. B 680, 34-43 (2009); CMS Collabora- 
tion, Phys. Lett. B 704, 411-433 (2011); ATLAS Collab- 
oration, Eur. Phys. J. C 71, 1647 (2011). 

[8] A. Falkowski, J. T. Ruderman, T. Volansky and J. Zu- 
pan, J. High Energy Phys. 05, 077 (2010); A. Falkowski, 
J. T. Ruderman, T. Volansky and J. Zupan, Phys. Rev. 
Lett. 105, 241801 (2010). 

[9] V. M. Abazov et al. (DO CoUaboration), Phys. Rev. Lett. 
105, 211802 (2010). 
[10] CMS CoUaboration, J. High Energy Phys. 07, 098 
(2011). 

[11] We use the convention that "momentum" refers to pc and 
"mass" to mc'^. 



[12] The CDF II detector uses a cylindrical coordinate system 
in which (p is the azimuthal angle, r is the radius from 
the nominal beam line, and z points in the proton beam 
direction. The transverse {r — <j)) plane is perpendicular to 
the z axis. Transverse momentum and energy are the re- 
spective projections of momentum measured in the track- 
ing system and energy measured in the calorimeter sys- 
tem onto the r ~ (j) plane, and are defined aa pr = p sin 6 
and Et = EsinO. Here, 9 is the polar angle measured 
with respect to the interaction vertex. Missing Et (JPt) 
is defined by ]^t = — JZi EtTh, where i is the calorimeter 
tower number for \ini\ < 3.6, and rii is a unit vector per- 
pendicular to the beam axis and pointing at the i^^ tower. 
The transverse mass rriT is defined as the invariant mass 
of the sum of the transverse projection of the momentum 
vector of one decay product and the missing energy vec- 
tor ^T- The pseudorapidity rj is defined as — ln(tan(^/2)), 
where 6 is measured with respect to the origin of the de- 
tector. We define the magnitude ^t = |^t|. 

[13] A. Abulencia et al., J. Phys. G 34, 2457 (2007). 

[14] F. Abe et al., Nucl. lustrum. Methods A 271, 387 (1988); 
D. Acosta et al. (CDF Collaboration), Phys. Rev. D 
71, 052003 (2005); The CDF-II Detector Technical De- 
sign Report, Fermilab-Pub-96/390-E; A. Abulencia et 
al. (CDF Collaboration), J. Phys. G 34, 2457 (2007); 
T. Aahonen et al. (CDF Collaboration), Phys. Rev. D 
77, 112001 (2008). 

[15] A. Sill et al., Nucl. lustrum. Methods A 447, 1 (2000); 
A. Affolder et al, Nucl. lustrum. Methods A 453, 84 
(2000); C.S. Hill et al, Nucl. lustrum. Methods A 530, 
1 (2000). 

[16] A. Affolder et al, Nucl. lustrum. Methods A 526, 249 
(2004). 

[17] L. Balka et al. Nucl. lustrum. Methods A 267, 272 



(1988). 

[18] S. Bertolucci et al. Nucl. Instrum. Methods A 267, 301 
(1988). 

[19] S. Kuhlmann et al, Nucl. Instrum. Methods A 518, 39, 
(2004). 

[20] G. Ascoh et al, Nud. Instrum. Methods A 268, 33 
(1988); A. Artikov et al, Nud. Instrum. Methods A 538, 
358 (2005). 

[21] D. Acosta et al (CDF CoUaboration), NucL Instrum. 
Methods A 461, 540 (2001); D. Acosta et al, Nucl. In- 
strum. Methods A 494 57 (2002). 

[22] D. Acosta et al (CDF Collaboration), Phys. Rev. Lett. 
94, 091803 (2005). 

[23] F. Abe et al (CDF Collaboration), Nucl. Instrum. Meth- 
ods A 271, 387 (1988). 

[24] S. Wilbur, "Search for Anomalous Production of Events 
with a. W or Z Boson and Additional Leptons", Ph.D. 
Thesis, University of Chicago (2011). 

[25] S. Paramonov, "A Limit on the Branching Ratio of the 
Flavor- Changing Top Quark Decay t -5> Zc." , Ph.D. The- 
sis, University of Chicago (2009). 

[26] T. Aaltonen et al (CDF Collaboration), arXiv:1004.1181 
[hep-ex] . 

[27] A. Abulencia et al (CDF CoUaboration), J. Phys. G 

Nucl. Part. Phys. 34 (2007). 
[28] D. Acosta et al (CDF Collaboration), Phys. Rev. D 72, 

032002 (2005). 

[29] M. L. Mangano, M. Moretti, F. Piccinini, R. Pittau, 
A. Polosa, "ALPGEN, a generator for hard multiparton 
processes in hadronic collisions", J. High Energy Phys. 
07, 001 (2003). 

[30] T. Sjostrand, S. Mrenna, P. Skands, "PYTHIA 6.4 
Physics and Manual", J. High Energy Phys. 05, 026 
(2006). 

[31] E. Gerchtein and M. Paulini, "CDF detector simulation 
framework and performance", arXiv:0306.031 [comp-ph]. 

[32] S. Hoeche et al., "Matching Parton Showers and Matrix 
Elements", arXiv:0602.031 [hep-ph]. 

[33] B. Cooper, Ph.D. thesis. University Coll. London (2006), 
FERMILAB-THESIS-2006-61. 

[34] M. Reece, private communication. 

[35] T. Junk, Nucl. Instrum. Methods A 434, 435 (1999). 



